perm filename REFLEC.F4[JC,MUS] blob
sn#078104 filedate 1973-12-17 generic text, type T, neo UTF8
00100 SUBROUTINE REFLECT(MAX,K)
00200 COMMON FREQ(3,0/50,50),FUNC(50),AMP(50),II(1),IJJ(3000)
00300 DO 100 N=0,MAX*2
00400 100 IF(FREQ(3,N,K).LT.0.0)FREQ(2,N,K)=-FREQ(2,N,K)
00500 DO 201 J=0,(MAX*2)-1
00600 DO 201 I=J+1,MAX*2
00700 IF(FREQ(1,J,K).OR.FREQ(1,I,K).EQ.99999.)GO TO 201
00800 IF(ABS(FREQ(1,J,K)).NE.ABS(FREQ(1,I,K)))GO TO 201
00900 IF(FREQ(1,J,K).GT.FREQ(1,I,K))GO TO 20
01000 C IF(ABS(FREQ(2,J,K)).LT.0.0001)FREQ(2,J,K)=0.0
01100 FREQ(2,I,K)=FREQ(2,I,K)-FREQ(2,J,K)
01200 FREQ(1,J,K)=99999.
01400 GO TO 201
01500 C20 IF(ABS(FREQ(2,I,K)).LT.0.0001)FREQ(2,I,K)=0.0
01600 20 FREQ(2,J,K)=FREQ(2,J,K)-FREQ(2,I,K)
01700 FREQ(1,I,K)=99999.
01800 201 CONTINUE
01900 RETURN
02000 END